Li 5 as a pseudorotating planar cluster 
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A pseudorotating state of Lis observed in a recent EPR experiment is investigated using the local 
density functional method. The calculated isotropic spin population indicates that a planar Civ 
structure is more consistent with the experimental result than a suggested trigonal bipyramid. A 
new pseudorotating state based on the planar structure is proposed. 

PACS numbers: 36.40. +d, 31.20.-d, 35.20.Bm, 35.20.Jv 



Despite a decade of theoretical investigation, the ge- 
ometric structures of small metallic clusters are poorly 
understood. This is partly due to the lack of a re- 
liable experimental technique and partly to the inade- 
quate accuracy of ab initio calculations. Only a few ge- 
ometries have been experimentally determined. Spectro- 
scopic data for the systems larger than the trimer are 
too complicated to extract structural information. On 
the other hand, highly accurate configuration-interaction 
(CI) calculations are too expensive even for small clus- 
ters when geometric optimization is required. A solu- 
tion to this challenging problem is to calculate spectro- 
scopic data or other observables using an accurate ab 
initio method for candidate geometries and to eliminate 
structures that are inconsistent with the experimental 
data. Currently available ab initio methods are accurate 
enough to distinguish quantitative differences in spectra 
due to different structures 0]. Recently, geometries of a 
few small lithium clusters have been partially determined 
in this way . 

For Li 3 , both matrix-isolated [j|]6) and cluster-beam |?|] 
electron paramagnetic resonance (EPR) spectra showed 
that all three nuclei are magnetically equivalent, indicat- 
ing that Li3 is an equilateral triangle [D^h]- On the other 
hand, theoretical calculations predicted two nearly de- 
generate energy minima, 2 Ax and 2 B 2 states, correspond- 
ing to C<i v isosceles triangles. This disagreement was 
explained by pseudorotation around the D^h symmetry 
point H||. Recently, extremum structures of the 2 A\ 
and 2 £>2 states were determined by high-resolution op- 
tical absorption spectroscopy fiofl . Furthermore, a frac- 
tional quantum number of pseudorotational angular mo- 
mentum due to Berry's phase 11 1^] has been observed 
in the rotational energy spectrum of Na 3 [ p^l4] ], sup- 
porting the pseudorotation ground state. 

A number of theoretical calculations p5|-^9[ predicted 
that lithium clusters are planar up to Lig. However, the 
optical absorption spectrum for Lig appears to be consis- 
tent with the theoretical spectra for a three-dimensional 
structure, eliminating the possibility of planar geometries 
jl],[| . Detailed geometrical parameters are still unknown. 
For Li7, EPR spectra ]2(J were interpreted as a pentag- 
onal bipyramid (D§h) in agreement with theoretical pre- 



diction (l^jr^,|l^,|T[ . As far as we know, no attempt to 
identify the structure of Lis has been made. 

Recently, Howard et al. [p2| investigated Lis m an 
adamantane matrix using EPR. The data indicate that 
all five nuclei are magnetically identical and this is in- 
terpreted as a pseudorotating cluster. Furthermore, the 
authors suggested that the ground state of Lis is not 
a planar but a three-dimensional structure, contrary to 
theoretical predictions. A trigonal bipyramid, which is 
nearly as stable as the planar structure 0,^1 , is a prime 
candidate for the most stable three-dimensional struc- 
ture. Since the pseudorotation mechanism for a trigonal 
bipyramid molecule such as PF5 is known ]24| , their spec- 
ulation seems reasonable. However, no conclusion can be 
made from the observed isotropic spin population (ISP). 
In order to extract the structural information out of the 
experimental data, quantitative theoretical calculations 
are desired. 

In this paper, we report a pseudopotential local den- 
sity functional calculation of Lis. Two low-energy iso- 
mers, a trigonal bipyramid ( 2 -Bi) and a planar structure 
( 2 Ai), are studied as candidates for the ground state. 
The isotropic spin population is computed in order to 
identify the structure observed by EPR. Furthermore, 
the possibility of two pseudorotation mechanisms is in- 
vestigated by searching the lowest potential barriers for 
permutation of the nuclei. 

The computational method used in this work is the 
same as the one previously used for other systems. The 
details are given in Ref. p5j ). The electron density is 
calculated using the local density approximation (LDA) 
[26 1 with a generalized norm-conserving pseudopotential 
[ 27 j . In order to increase transferability, a fractionally oc- 
cupied electronic configuration, ls 2 2s 08 2p - 2 , is used to 
construct a pseudopotential ]27y . A core radius R c = 2.1 
is sufficiently accurate for all orbitals ^3|. Wave func- 
tions are expanded in a plane wave basis set using a sim- 
ple cubic supercell. Since Li 3 has a relatively large dipole 
moment, a large lattice constant, a = 60 a.u., was used 
to avoid a long-range dipole interaction between repeated 
cells. For other sizes, a — 30 a.u. was sufficiently large. 
The same cutoff energy E cut = 11.2 Ry was used for all 
sizes. This requires more than 130000 plane waves for the 
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largest cell. Both electronic and geometric structures are 
optimized by the steepest descent method. The ISP at 
the zth nuclear position, Ri is evaluated as 
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where V>n,<r is the n-th molecular orbital with a spin a = 
±1 and ip2s is the 2s orbital of an isolated Li atom. All 
calculations have been carried out using the CM-200 and 
CM-5 computers. 

We have calculated Li 3 , Li 5 and Li7 clusters at the 
same level of accuracy. The ISP's for Li 3 and Li7 is com- 
pared with the experimental data in order to provide the 
degree of accuracy of our calculation for Lis ■ 

For Li3, two equilibrium triangle structures, 2 B 2 (an 
obtuse isosceles triangle) and 2 A\ (an acute isosceles tri- 
angle) are nearly degenerate as predicted in other calcu- 
lations Hll. The 2 B 2 state appeared to be more stable 
than the A\ state. However, the difference in total en- 
ergy is only 8.3 meV. The experimental value, 3.2 meV 
JlO| |, is within the accuracy of the present calculation. 
The conical intersection of two Born-Oppenheimer (BO) 
surfaces at the D^h symmetric point is only 66 meV (the 
experimental value = 33 meV) above the 2 B 2 minimum, 
indicating that the system is quantum mechanically de- 
localized and forms a D%h instead of a C 2v structure. 
These results are consistent with the pseudorotating state 
identified by the EPR spectroscopy The ISP's for 

both extremum geometries are listed in Table |. The 
total ISP, ptotai is 0.48 and 0.78 for the 2 A X and 2 B 2 
states, respectively. Since the pseudorotating state oc- 
cupies both extremum states, the actual ptotai should be 
between these values. Assuming equal weights on the two 
states, ptotai—^- 63 is obtained, in good agreement with 
the experimental value, p totQ /=0.69. 

The EPR experiment determines the number of mag- 
netically equivalent nuclei and the isotropic spin pop- 
ulation at each nucleus but it is not capable of deter- 
mining the detailed geometry such as equilibrium bond 
distances. Recently, the extremum geometries were pre- 
cisely determined using high-resolution photoabsorption 
spectroscopy through the analysis of vibronic spectra 
p0| . The present calculation gives excellent agreement 
with the experiment, as shown in Table H. 

Beyond the trimers, both optical and EPR spectra are 
so complicated that the assignment or interpretation of 
spectra is almost impossible without theoretical assis- 
tance. For L17, however, two distinct groups of nuclei, 
two nuclei with large ISP (pi=0.25) and five nuclei with 
very small ISP (p 2 =-0.015) were found from EPR |o), 
consistent with the theoretically predicted pentagonal 
bipyramid [ p^|p^Jl9| , pl[ | . The present calculation predicts 
that two apical nuclei have pi=0.23 and five nuclei form- 
ing a pentagonal ring carry p 2 —— 0.01, in good agreement 
with the experiment. 

The results for Li 3 and Li7 suggest that these calcula- 
tions are sufficiently accurate to predict the ISP of small 



lithium clusters. Therefore, the geometry of small clus- 
ters can be determined by comparing theoretical ISP's 
with experimental data. 

For Lis , most theoretical calculations Jl5|-|l9|l predicted 
that the 2 A\ state of a C 2v planar structure [See Fig. 0(a)] 
is the ground state. However, in contradiction to the 
previous calculations, we found that a distorted trigonal 
bipyramid with the C 2v symmetry [See Fig. ^(a)] is lower 
in total energy than the 2 Ai state within the LDA. This 
result seems consistent with the speculation by Howard 
et al. |22|. However, the energy difference is only 46 
meV/atom. Since the accuracy in LDA total energy is 
not always satisfactory |28| , determining the ground state 
by the total energy alone is not reliable if there are nearly 
degenerate states. 

In Table III, the theoretical ISP values of two struc- 



tures are compared with the experimental value. The to- 
tal ISP, ptotai = 0.53, for the trigonal bipyramid is much 
smaller than the experimental value, 0.71. On the other 
hand, ptotai — 0.72 for the planar structure nearly coin- 
cides with the experimental value. This implies that the 
observed Lis cluster is the planar cluster. The small ISP 
for the trigonal bipyramid is due to the large population 
of the p-character electron which stabilizes this geometry 
significantly. Large sp mixing is particular to Li clusters 
and does not appear in other alkali-metal atom clusters. 

The experiment also showed that the five nuclei are 
magnetically equivalent. Only a pentagonal ring (D^h) 
has five equivalent nuclei without dynamical transforma- 
tion. However, it is an unstable state with total energy 
of 1.06 eV above the the lowest 2 B\ state and, therefore, 
energetically not accessible. The observed spectra must 
be due to the rapid transformation of one of the low en- 
ergy structures. In this paper, we assume the system is 
pseudorotating by quantum-mechanical tunneling. 

By analogy from the Berry pseudorotation of PF5, the 
trigonal bipyramid permutes its nuclei through the C± v 
transition state, as illustrated in Fig. ^|. For Lis, however, 
since the unpaired electron has a node on the atoms 3-1- 
3 in Fig. ||(a) and on atoms 2-1-2 in Fig. ||](c), the elec- 
tronic wave functions before and after the transition are 
orthogonal. Therefore, tunneling through the C± v point 
is essentially prohibited. Furthermore, the potential bar- 
rier height, 0.33 eV is too high to allow a rapid tunneling 
at low temperature. The pseudorotation of the trigonal 
bipyramid is not likely to occur within the experimental 
time scale. 

Similarly to the pseudorotation of Li 3 , we expect that 
the planar Lis pseudorotates around the D$h symmetric 
point if the potential barrier to the cyclic permutation 
of nuclei is sufficiently low. The D$h point is 0.82 eV 
above the 2 Ai state, which is too high to pass through 
even if the zero-point energy is taken into account. How- 
ever, we found a relatively low transition state to reach 
a nearest equivalent extremum in the planar structure. 
The path from one extremum to another is illustrated in 
Fig. |l|. The barrier height, 0.18 eV, is higher than that 
of Li 3 but significantly lower than that of the trigonal 
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bipyramid pseudorotation path. In contrast to the trigo- 
nal bipyramid, the nodal planes before and after tunnel- 
ing are nearly parallel and thus the tunneling is allowed. 
Although a nearly free pseudorotation like the Li3 case is 
not possible for Lis, the tunneling probability of the pla- 
nar structure is expected to be much larger than that of 
the trigonal bipyramid, supporting the idea that the pla- 
nar structure is the ground state and in consistent with 
the CI calculation (Table |v|). 

In order to make all nuclei equivalent, tunneling has 
to be faster than the observation time of the experiment. 
Unfortunately, a quantum Monte Carlo simulation is nec- 
essary to determine the quantum-mechanical behavior of 
Li nuclei. A recent quantum Monte Carlo simulation 
p9| indicates that zero-point energy and tunneling play 
a crucial role even for large clusters. Furthermore, the es- 
timated zero-point energy for the bulk Li crystal is as big 
as 33 meV/atom |3(J. Since lithium is the third lightest 
element, the quantum-mechanical behavior of the nuclei 
is not surprising. 

The dimension of the pseudorotation is 6 for Li 5 in 
contrast to 2 for Li 3 . Because of this high dimensional- 
ity, the Lis pseudorotation is not straightforward. The 
permutation path described in Fig. |l] is actually not a 
simple rotation. In terms of rotational motion, the tran- 
sition is not between the nearest neighbors as shown in 
Fig. U due to the high dimensionality of the system. The 
path is a closed loop in the six-dimensional space and 
its projection onto a two-dimensional plane rotates three 
times (67r) around the symmetry point. It is inter- 
esting to calculate Berry's phase associated with the ge- 
ometrical transformation along the path |L2 31 and the 
fractional quantization of rotational motion. However, 
calculation of Berry's phase is not simple for Lis because 
of its high dimensionality and the degeneracy at the D^h 
point @|l|- 

In conclusion, the theoretical isotropic spin population 
indicates that the planar C'2v{ 2 ^i) structure is more con- 
sistent with the EPR data than the C2 V ( 2 Bi) trigonal 
bipyramid. Furthermore, the pseudorotation of the trig- 
onal bipyramid is expected to be very slow due to the 
symmetry of the electronic wave function and the high 
potential barrier. A new pseudorotation path thorough 
a transition state at C2 V ( 2 Bi) is proposed for the planar 
structure. The low potential barrier height probably per- 
mits rapid tunneling between extrema in the BO surface. 
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TABLE I. Isotropic spin populations for Li3. The experi- 
mental data are for the pseudorotating state. 
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a Reference 6 Reference || c This work 
d References [§-[7| 

TABLE II. Comparison of calculated extremal geometries 
of Li3. The cohesive energy E co h is in eV, relative energy to 
the 2 £?2 state AE in eV, side length R in a.u., and apex angle 
6 in degrees. 
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TABLE III. Isotropic spin populations p; for Lis. See Figs. 
[0(a) and @(a) for the location of the ith nucleus. 
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TABLE IV. 


Comparison 


of theoretical calculations for the 


low-energy geometries and cohesive energ 
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FIG. 1. Pseudorotation path for the planar Lis ( 2 Ai). The 
center structure corresponds to the transition state between 
two extremum states (the left and the right figures). 
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FIG. 2. Pseudorotation path for the trigonal bipyramid 
( 2 Bi). The center structure corresponds to the transition 
state between two extremum states (the left and the right 
figures) 




FIG. 3. Pseudorotation around the D^h symmetry point. 
Arrows represent the transition shown in Fig. hi 
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